Hyper-domains in exchange bias micro-stripe pattern 
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A combination of experimental techniques, e.g. vector-MOKE magnetometry, Kerr microscopy 
and polarized neutron reflectometry, was applied to study the field induced evolution of the magneti- 
zation distribution over a periodic pattern of alternating exchange bias stripes. The lateral structure 
is imprinted into a continuous ferromagnetic/antiferromagnetic exchange-bias bi-layer via laterally 
selective exposure to He-ion irradiation in an applied field. This creates an alternating frozen-in in- 
terfacial exchange bias field competing with the external field in the course of the re-magnetization. 
It was found that in a magnetic field applied at an angle with respect to the exchange bias axis 
parallel to the stripes the re-magnetization process proceeds via a variety of different stages. They 
include coherent rotation of magnetization towards the exchange bias axis, precipitation of small 
random (ripple) domains, formation of a stripe-like alternation of the magnetization, and devel- 
opment of a state in which the magnetization forms large hyper-domains comprising a number of 
stripes. Each of those magnetic states is quantitatively characterized via the comprehensive analysis 
of data on specular and off-specular polarized neutron reflectivity. The results are discussed within 
a phenomenological model containing a few parameters which can readily be controlled by designing 
systems with a desired configuration of magnetic moments of micro- and nano-elements. 

PACS numbers: 75.60.Ch; 75.60.Ej; 75.75.+a 61. 12. Ha; 



I. INTRODUCTION 



The exchange bias (EB) effect, which is due to ex- 
change coupling between ferromagnetic (F) and antifer- 
romagnetic (AF) layers, is expressed via a shifted hys- 
teresis loop away from zero field. The shift is attributed 
to the frozen-in global unidirectional anisotropy of the 
system. Due to its intriguing physics and importance for 
device application the EB effect is persistently under ex- 
tensive study (see Refsji^^). Spacial alteration of the 
EB field brings qualitatively new physics into EB sys- 
tems and creates new promising technological perspec- 
tives. Deep understanding of the role of competing in- 
teractions in this class of materials is required e.g. to 
manufacture F/AF films with imprinted in-plane ferro- 
magnetic domain pattern with desired morphology. One 
of the motivations is to design planar magnetic pattern- 
ing of a continuous film as an alternative to magnetic 



grains, clusters or structures patterned by lithographic 
methods. With such systems one could avoid the problem 
of very small feature sizes, where the long-term thermal 
stability of the magnetic elements is lost due to super- 
paramagnetic fluctuations^. 

The EB effect is usually set via cooling the system 
below the blocking temperature of the AF layer in a 
magnetic field saturating the ferromagnetic counterpart. 
The size and the sign of the EB effect depends on the 
choice of F/AF materials in contact and often can be 
altered via changing the field cooling protocol 7 *^. On 
the other hand, the EB field direction and strength can 
be selectively altered by ion bombardment of the F/AF 
bi-layer subjected into an external field of the opposite 
directio n 10 ' 11 . Depositing a grid protecting some areas of 
the sample one can preserve the EB field in those areas, 
while altering its direction in the unprotected regions. 
The method of Ion Beam Induced Magnetic Patterning 
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(IBMP)A2, opens wide perspectives to produce various EB 
patterns. Here we report on the magnetic properties of 
a IBMP produced stripe-like pattern with the EB-field 
set antiparallel in neighboring stripes so that the net EB 
field is compensated and the system should posses global 
uniaxial, instead of unidirectional, anisotropy. Therefore 
it is expected that the ferromagnetic stripe domains in 
such a system display alternating magnetization direc- 
tions in the remanent demagnetized state of the sample. 
Details on the sample preparation and the results of ex- 
perimental and theoretical studies of the magnetization 
reversal mechanism of the system for the field applied 
along the exchange bias axis can be found in our recent 
publication 13 . There by use of the magneto-optical Kerr- 
effect (MOKE) in vector-MOKE configuration, Kerr mi- 
croscopy and polarized neutron reflectometry (PNR) it 
has been shown that the system exhibits a hysteresis rich 
in details and with a complex domain structure. Surpris- 
ingly, it was found that in the easy axis configuration the 
magnetic state after reversal of one of the both types of 
stripe domains never shows a clear antiparallel domain 
structure. Instead, at low fields the magnetization in the 
different EB stripes is periodically canted with respect to 
the EB axis so that alternating antiparallel ordering of 
domain magnetization projections onto the stripe axis co- 
exists with a macroscopic magnetic moment perpendicu- 
lar to the anisotropy axis. This effect is explained within 
the framework of a simple phcnomenological model which 
takes into account competing interfacial and intralayer 
exchange interactions. According to the model, within a 
certain range of parameters including, e.g. the ratio be- 
tween ferromagnetic layer thickness and the stripe width, 
the system reveals an instability with respect to the tilt of 
magnetization to the left, or to the right against the EB 
stripe induced anisotropy axis. In our previous studyi 3 - 
it is always found flipped only in one of two nominally 
equal directions so that the system always carried an ap- 
preciable residual magnetization not collinear with the 
anisotropy axis and the magnetic field. 

The finding in Ref^ 3 - may have more general and far 
going consequences for the understanding of the inher- 
ent physics of the EB effect. Areas with alternating EB 
fields must exist in any F/AF bi-layers coupled via ex- 
change interaction through a common atomically rough 
interfac o 14 ! 15 . Alternating interfacial fields in this case 
are generated by uncompensated spins in AF areas in 
neighboring atomic steps which are shifted up or down 
with respect to each other for half of the magnetic unit 
cell of the Ising type antiferromagnet. Those interfa- 
cial fields randomly alternate over the F/AF interface 
and compete with the lateral exchange interaction which 
favors a homogeneously magnetized ferromagnetic film. 
Depending on film thickness and strength of the inter- 
actions, e.g. the interfacial F/AF vs lateral exchange in 
the ferromagnetic film, the competition may result in a 
state with magnetization of the Heisenberg ferromagnet 
(inhomogeneously) tilted away from the external field ap- 
plied collinearly with EB direction. In view of that the 



magnetization distribution in the model system with a 
controlled EB field alternation deserves further compre- 
hensive consideration. 

First of all we admit that the previously observed^ 
preferential large tilt of the net magnetization away from 
the symmetry axis can be explained by a tiny misalign- 
ment between frozen-in EB fields in irradiated and pro- 
tected stripes. Such a misalignment determines a pref- 
erential direction of the small net EB field noncollinear 
with stripes. In the vicinity of the instability point the 
net magnetization naturally appears in the direction of 
this field. It is quite a challenging technological task to 
set both EB fields exactly antiparallel and collinear with 
the stripes. This is not a goal of the present paper, where 
instead, we report on the results avoiding this problem 
by an deliberate tilt of the external field direction by an 
angle as high as 45° regarding to the anisotropy axis set 
along the stripes. Then a little misalignment between the 
EB fields in neighboring stripes plays a minor role. Mea- 
surements in an asymmetric regime, on the other hand, 
disclose many details on domain organization and evo- 
lution which are otherwise hidden, but absolutely cru- 
cial for a complete understanding of the re-magnetization 
mechanisms in EB patterned arrays and other types of 
systems with alternating EB fields. The bulk of infor- 
mation is mostly gained via the quantitative analysis of 
data on PNR. The scattering intensity distribution was 
measured over a broad range of angles of incidence and 
scattering and accomplished with a full polarization anal- 
ysis at different magnetic fields along the hysteresis loop. 
For fitting the specular reflectivity data we used an orig- 
inally developed least-squares software package^ which 
allows for simultaneous evaluation of all four measured 
reflectivities in one cycle. We compare the results of 
our fits to vector-MOKE measurements. For supporting 
the interpretation of our data on off-specular diffuse and 
Bragg scattering Kerr-microscopy (KM) images^ were 
also taken. 



II. EXPERIMENTAL DETAILS 

The sample studied is a TaO-Ta(8.7 nm)/ 
Co 7 oFe 30 (28.0 nm)/ Mn 83 Iri 7 (15.2 nm)/ Cu(28.4 nm)/ 
SiO 2 (50.5 nm)/ Si(lll) film stack. The initial EB 
direction was set by field cooling in a magnetic in-plane 
field of 100 mT after an annealing step for 1 h at 
275 °C which is above the blocking temperature of the 
antiferromagnetic material. Subsequently, IBMP using 
Hc + ions with a fluency of 1 x 10 15 ions/cm 2 at 10 
keV was applied at a magnetic field of 100 mT aligned 
opposite to the initial EB direction. This resulted in a 
stripe-like domain pattern of equally spaced stripes with 
a width of 2.5 /im and a periodicity of A = 5 /im and 
with alternate sign of the unidirectional anisotropy, and 
hence of the EB in neighboring stripes^ 

The evolution of magnetization as a function of ap- 
plied field was recorded with a vector-MOKE setup as 
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described in Ref. 18 . Both magnetization projections were 
accessed via measuring the Kerr angle parallel and per- 
pendicular to the field. The projections were measured 
in the longitudinal MOKE configuration, i.e. with the in- 
cident light linearly polarized within the reflection plane. 
For determination of the magnetization component per- 
pendicular to the field the sample and the magnetic field 
direction were rotated simultaneously by 90° about the 
normal to the surface. Assuming that the Kerr an- 
gles, 8^ oc Mcos 7 /, measuring the magnetization vec- 
tor projection and the Kerr angle parallel to the field, 
and oc Msm.7, corresponding to the perpendicular 
magnetization component, have the same proportionality 
coefficients, one can determine the tilt angle 7 through 
tan 7 = #k/^k an d the normalized length of the mag- 
netization vector M/M a = \J {O 1 ^) 2 + (O^f/O^ with M s 

being the modulus of the magnetization and 9 S K the Kerr 
angle, both in saturation. This allows to completely de- 
termine the direction of the mean magnetization vector 
M and its absolute value M = \M\ reduced due to do- 
mains. In different domains the magnetization vector 
deviates by angles A7 from the direction of the mean 
magnetization and hence is tilted by 7 = 7 + A7 against 
the applied field. The mean angle 7 is determined by 
the equation (sinA7) co h = 0, where the averaging runs 
over the spot coherently illuminated by the laser beam. 
Then the normalized longitudinal and, correspondingly, 
transverse MOKE signals can be written as: 

(cos7)coh = M/M s cos 7 = (cos A7) coh cos7, (1) 
(sin7) coh = M/M s sin 7 = (cos A7) coh sin 7, (2) 

where the mean magnetization M = M s (cos A7) co h. 

Further insight into the microscopic rearrangement 
of magnetization was achieved by a high resolution 
magneto-optical Kerr microscope (KM}i£ that is sensi- 
tive to directions orthogonal to the fielcU^. 

Neutron scattering experiments were carried out with 
the HADAS reflectometer at the FRJ-2 reactor in Jiilich, 
Germany. Details of the measuring geometry can be 
found in Ref^. In the present experiment the sample 
was aligned so that the field orientation was tilted by 45° 
against the EB axis. The monochromatic neutron beam 
with the wavelength of 0.452 nm incident onto the sam- 
ple surface under the glancing angle a>i is scattered under 
the angle a/, so that for specular reflection 07 = ctj. The 
incident polarization vector Pi was directed either paral- 
lel or antiparallel to the magnetic field and perpendicular 
to the scattering plane. In each of the two directions of 
Pi the final spin state was analyzed with respect to the 
same direction with an efficiency Pf — \Pf\ < 1. 

Specular PNR provides information similar but not 
identical to vector-MOKE magnetometry. Both meth- 
ods probe projections of the magnetization vector aver- 
aged over its deviations due to magnetic domains and 
other inhomogeneities within the coherence volume of 
photons, or neutrons, respectively. In the case of MOKE 
the laser beam is rather coherent all over the isotropic 



light spot illuminating the sample surface. This is not 
the case for PNR. The neutron source is essentially in- 
coherent, but the neutron beam is well collimated in the 
reflection plane, while the collimation is usually relaxed 
perpendicular to this plane. Hence the cross section of 
the coherence volume with the reflecting surface is rep- 
resented by an ellipsoid with dramatically extended axis 
along the beam projection onto the surface. At shallow 
incidence this extension, i.e. the longitudinal coherence 
length, may spread up to some fractions of a millime- 
ter. In contrast, the other ellipsoid axis, i.e. the coher- 
ence length across the incoming beam is short and only 
amounts to about 10 nm. Therefore, the coherence 2D 
ellipsoid covers only a very small area of the sample, and 
the observed PNR signal is a result of two sorts over aver- 
aging. The first one includes a coherent averaging of the 
reflection potential, e.g. over directions of the magneti- 
zation in small (periodic and random) domains, and runs 
over the coherence area. Second, the reflected intensities 
from each of those areas are summed up incoherently. 

If the mean magnetization averaged over the coher- 
ence area is nonzero and collinear with the external 
field (applied parallel to the neutron polarization axis 
and perpendicular to the scattering plane) then spec- 
ular reflection does not alter the neutron spin states 
and only two non-spin-flip (NSF) reflection coefficients 
1Z ++ 7^ 1Z are measured, while both spin-flip (SF) 
reflectivities K + ~ = TZ~ + = 0. In this case NSF re- 
flectivities are uniquely determined by the mean optical 
potential, e.g. by the mean projection of the magneti- 
zation proportional to (cosA7) co h, where A7 is the tilt 
angle of the magnetization in domains smaller than the 
coherence length 2 ^. 

If the mean magnetization direction makes an angle 
7 = (7 — A7) co h with the applied field then the SF re- 
flectivities 1Z + ~ = TIT + ^ and are proportional to 
sin 2 7, i.e. to the mean square of the magnetization pro- 
jection normal to the field. At the same time the dif- 
ference, 1Z + + — 1Z , between the NSF reflectivities is 
proportional to cos 7, i.e. to the projection of the mean 
magnetization within the coherence volume onto the ap- 
plied field direction. Crossing a number of small periodic 
and random domains in only one direction the coher- 
ence ellipsoid still covers a very small area of the sample. 
Therefore, the angle 7 may vary along the sample surface 
and reflectivities have to be incoherently averaged over 
the spread of 7. It is important to note that NSF and 
SF reflectivities are complicated nonlinear functions of 
(cosA7) co h, which may also vary along the sample sur- 
face. Hence, such an averaging is, in general, not a trivial 
procedure. 

If the value of (cos A7) co h is, however, unique^ for the 
whole sample surface, then 

TL++ - TT- oc (cos7) inc (3) 
Tl+- = TZ-+ oc (sin 2 7) inc (4) 

are respectively proportional to cos 7 and sin 7 inco- 
herently averaged over the reflecting surface with the 
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proportionality coefficients nonlinearly depending on 
(cosA7) co h- Because of nonlinearity the parameters 
(cos7)inc, (sin 2 7)j nc and (cosA7) co h can only be found 
via fitting of the data for all NSF and SF reflectivities. 
After that one can determine the mean value (cos 7) w 
(cos7~)i llc (cos A7) co h, under the condition: (sinA7) co h = 
0. This constraint is similar to that applied above for 
vector-MOKE and hence the mean magnetization projec- 
tion onto the field direction is in this case equally mea- 
sured by both methods: MOKE and PNR. Then agree- 
ment between results of PNR and vector-MOKE can be 
used to prove the assumption above. Alternatively, PNR 
is able to deliver an important information complement- 
ing vector-MOKE data. 

Due to the strong anisotropy of the coherence ellip- 
soid PNR can probe a variation of 7 in the correspond- 
ing direction measuring fluctuations of the magnetization 
not accessible for MOKE. In particular, with PNR one 
gains a direct access to the dispersion E 2 = (sin 2 7)i nc — 
(sin7~) 2 nc > which quantifies those fluctuations. If, for 
instance, E = these fluctuations are absent, then at 
(sin 2 7)i nc = (sin7~) 2 nc ^ the homogeneous sample mag- 
netization is homogeneously tilted by the angle 7 to the 
left, or to the right with respect to the applied field. PNR 
is not sensitive to the sign of 7, which can be determined 
via vector-MOKE. On the other hand, MOKE cannot 
provide any information about, e.g. the totally demag- 
netized structure when (sin 7} = (cos 7} = 0. In this 
case missing information can readily be retrieved from 
the PNR data. This can already be seen considering two 
limiting situations. The limiting value E = 1 is achieved 
in the state with large domains where the magnetization 
is perpendicular to the applied field. The other limit 
E = is reached if demagnetization occurs on a scale 
smaller than the coherence area. In the latter case spec- 
ular reflection is accompanied by off-specular scattering. 

Off-specular PNR can, in contrast to MOKE, directly 
measure the spread of magnetization directions due to 
domains crossed with the long axis of the coherence el- 
lipsoid. Periodic deviations A7 give rise to Bragg diffrac- 
tion, while random fluctuations cause diffuse scattering. 
The positions of the Bragg peaks in the reciprocal space 
determine the period of the domain structure along the 
largest coherence axis, while the extension of diffuse scat- 
tering is due to the correlations of the random compo- 
nent of the magnetization. Via fitting of intensities of off- 
specular scattering one can deduce the magnetization dis- 
tribution between neighboring stripe domains and in ad- 
dition to (cos A7) co h to determine one more parameter, 
(sin 2 A7) co h, characterizing magnetization fluctuations 
due to random ripple domains. Then one can also calcu- 
late the dispersion a 2 — (cos 2 A7) co h — ( cos ^7)coh — 
quantifying the microscopic arrangement of magnetiza- 
tion fluctuations within the coherence length. 

Often, the set of parameters indicated above is not suf- 
ficient to describe experimental data of PNR and to infer 
from it a domain state. Such a situation takes place when 
the mean value (cos A7) co h averaged over the coherence 



range varies along the sample surface. Then incoherent 
averaging must take into account that specular reflection 
and off-specular scattering occur from areas with differ- 
ent optical potentials. As we shall see below, within a 
certain range of applied fields the magnetization in our 
system evolves via formation of large (hyper-) domains 
comprising a number of small stripe and ripple domains. 
In contrast to the case of the conventional domain state, 
now the magnetization in different hyper-domains differs 
not only in directions, but also in absolute values. This 
is due to the fact that each type of hyper-domains is 
characterized by a particular arrangement of the mag- 
netization over stripe and ripple domains belonging to 
this type. Each of them acquires its own set of pa- 
rameters, e.g. (cos A7y} co h, (sin 2 A7j) coh , (cos^)^ and 
(sin 2 7 ■) inc , where the superscript j indexes the type of 
hyper-domain. Hf^, and SF, Ti-f T , reflectivities, as well 
as cross sections of off-specular scattering, are also to be 
indexed accordingly. Then the weights Wj of different 
types of hyper-domains can be determined along with 
the above listed parameters via the fitting of experimen- 
tal data. This would allow to totally characterize the 
magnetic states of the system passing through along the 
hysteresis loop. 

III. EXPERIMENTAL RESULTS 

A. Vector-MOKE 

The evolution of two Cartesian projections of the nor- 
malized to saturation mean magnetization vector m = 
M/M sat determined with vector-MOKE is depicted in 
Fig. Ha) and (b). Fig. 12(c) and (d) illustrate the field 
dependence of the absolute value m = \m\ of the nor- 
malized magnetization vector m and its tilt angle 7 with 
respect the field applied at the angle of 45° relative to 
the stripes axis. Comparing these plots one can admit 
several stages of the re-magnetization process as visual- 
ized in Fig. [2 Reduction of the negative field from sat- 
uration leads firstly to relatively slow coherent rotation 
of the magnetization vector away from the field direc- 
tion while maintaining its absolute value [Fig. [2] (a)]. At 
small positive field H = H c \ « 0.6 mT this process is 
suddenly interrupted apparently due to domain forma- 
tion [Fig. [2] (b)] . This stage of the process is completed 
at H ps +1.4 mT when the magnetization loss is about 
half of its magnitude. At further increase of the applied 
field the magnetization partially restores its magnitude 
up to almost 63% of the nominal value. At the same 
time the mean magnetization vector is directed at an an- 
gle 7 » 35° with respect to the applied field [Fig. [2] (c)]. 
Within quite a broad interval of fields the system stays 
in a relatively rigid state with the mean magnetization 
directed almost normal to the stripes. The next dramatic 
event occurs at H c 2 between H w 5.5 mT and H w 7.0 
mT when the magnetization again looses and partially re- 
stores its absolute value [Fig. [21(d)]. Within this stage the 
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FIG. 1: (Color online) Vector-MOKE data and results of fits 
to the PNR data, (a) Normalized hysteresis loop of the longi- 
tudinal magnetization component, (b) normalized hysteresis 
loop of the transverse magnetization component, (c) normal- 
ized value of the length of the magnetization vector, and (d) 
mean value of the angle between magnetization and external 
field. The field is applied at 45° with respect to the EB axis. 
Lines represent MOKE and symbols PNR data. The symbols 
represent the results of the the weighted average between first 
and second domain (see text). The ascending curve is drawn 
in color and the descending one in black. 



vector m rapidly starts to rotate towards the magnetic 
field direction and is aligned along the field at H « 6.0 
mT. At higher fields the magnetization, surprisingly, con- 
tinues to rotate further away from the field direction and 
at H « 6.5 mT it arrives at a maximal tilt angle 7 ~ 9°. 
This stage of the the re-magnetization process is appar- 
ently governed by a domain re-arrangement mechanism 
which restores the magnetization absolute value up to 
about 65% of its nominal. Rapid restoration of magneti- 
zation follows, presumably, through an additional inter- 
mediate stage of domain evolution within the field frame 
7.0< H<H sat [Fig. 1(e)]. 

The descending branch of the hysteresis loop [the first 
stage is shown in Fig. [2] (f)] generally repeats all main 
steps of the magnetization evolution. However, the mag- 
netization vector does not passes them in the exactly 
reversed order, as would be expected. Instead, along the 
descending branch the vector M continues to rotate in 
the same direction passing by the state with the magne- 
tization along the field as in the case of the ascending 
branch of the hysteresis loop. Finally the vector m ac- 
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FIG. 2: (Color online) Evolution of the mean magnetization 
along the MOKE hysteresis loop. 



complishes the full circle of 360° and then slightly rotates 
back approaching negative saturation. The intrinsic rea- 
son of such a behavior should find its explanation below 
after more detailed analysis of the complete scope of the 
data. 

Here we just mention that the hysteresis loops are 
shifted exhibiting a global EB effect and hence a residual 
frozen-in magnetic field. This indicates an incomplete 
compensation of fields frozen-in different sets of stripes. 
Further insight into the arrangement of magnetization 
over stripes is gained by use of Kerr microscopy. 



B. Kerr microscopy 

A sequence of KM images taken along the hysteresis 
loop for the present sample with the field applied along 
the exchange bias axis, perpendicular, and at 45° to 
this axis were briefly reported recently^. It was demon- 
strated, that in the latter case the evolution of the mag- 
netization distribution recorded in the images, e.g. in 
those presented in Fig. [3l elucidates the role of various 
mechanisms involved in the re-magnetization process ac- 
cording to the typical stages of the re-magnetization pro- 
cess shown in Figs. [2] (a-f). Figs. [3ja-e) illustrate the re- 
magnetization scenario along the ascending branch after 
saturation in a negative field. Fig.[3]Ja) shows a measure- 
ment performed at a small positive field of 0.3 mT. At 
this field the magnetization already appreciably deviates 
from the field axis but it is not yet reversed in neither 
of the stripes. The periodic KM contrast indicates an 
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FIG. 3: Kerr microscopy images taken during the magnetiza- 
tion reversal with the field aligned diagonally to the EB axis. 
The images were recorded with a mixed Kerr sensitivity tuned 
nearly to the transverse magnetization direction. 



angle between the directions of magnetization in neigh- 
boring stripes. Some rippling, predominantly in the He + 
bombarded regions, can also be observed. At 0.6 mT 
[Fig. O^b)] the reversal occurs through the formation of 
large hyper-domains separated along one of the stripes. 
In one of such hyper-domains the magnetization of the 
stripes is not yet reversed. In the other hyper-domain the 
magnetization in one set of stripes, i.e. in this case in the 
He + bombarded, is reversed as seen via a large optical 
contrast. At 1.4 mT no hyper-domains are seen and the 
magnetization projections onto the EB axis in neighbor- 
ing stripes seems to be aligned antiparallel [Fig. [3] (c)]. 
Further increase of the applied field changes the scenario 
of the re-magnetization process, as seen in Figs. [3J (d) 
and (e). Now it proceeds through gradual decrease of 
the width of stripes with the most unfavorable magne- 
tization direction. Fig. [3J (f ) was taken in the backward 
branch. It shows that reappearance of the continuous 
stripes with a negative projection onto the field direc- 
tion is preceded by precipitation of small ripple domains. 
Further decrease of the positive field and its subsequent 
alternation restores the periodic structure via coalescence 
of ripple domains in corresponding sets of stripes^ 7 . This 
process, however, does not occur simultaneously all over 
the sample surface, but again involves the formation of 
hyper-domains. Some of them contain homogeneously 
magnetized stripes, while in the others the magnetiza- 
tion of one set of stripes is broken into ripple domains. 



It should be noted that KM images sample, in contrast 
to vector-MOKE and PNR, quite a restricted area of 
the surface. They provide, however, a rather solid back- 
ground necessary for a quantitative analysis of PNR data 
and consequent characterization of magnetization distri- 
bution between stripe, ripple and hyper-domains over the 
whole sample. 
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FIG. 4: (Color online) Polarized neutron reflectivity mea- 
surements performed at different magnetic fields applied at 
the angle x — 45° with respect to the EB axis. The symbols 
present measurements of non-spin flip reflectivities 1Z + + and 
7Z and spin flip reflectivities 1Z + ~ and TV + . The lines 
represent fits to the data points. 



C. Neutron scattering 

Neutron scattering experiments were performed using 
a position sensitive detector (PSD). The PSD records, ad- 
ditionally to the specular reflection from the mean neu- 
tron optical potential, magnetic Bragg diffraction from 
the periodic stripe array and off-specular diffuse scat- 
tering from domains smaller than the long axis of the 
coherence ellipsoid. The data were taken at several posi- 
tive field values of the hysteresis. The magnetic field was 
kept parallel to the field guiding the neutron polariza- 
tion in order to avoid neutron depolarization. Prior to 
the measurements, the sample was saturated in a nega- 
tive field. Specular reflectivities were extracted from the 
PSD maps. In Fig.[2J several representative experimental 
curves together with fits to the data are displayed. Most 
of the presented data are collected at fields correspond- 
ing to the KM images in Fig. [3j This allows for a qual- 
itative interpretation of the specular reflectivity curves 
along with the KM measurements in Fig. [3J 

At 0.3 mT splitting of the NSF reflectivities with TZ~ ~ 
being higher than 1Z + + and considerable SF reflectivities 
indicate an appreciable tilt of the net magnetization, al- 
most homogeneous in accordance with Fig. [3J (a), away 
from the direction antiparallel to the applied field. 

At 0.8 mT the SF reflectivity is slightly increased and 
the splitting of the NSF reflectivities is reduced compared 
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to 0.3 mT. This is attributed to a further increasing tilt 
of magnetization and reduction of its absolute value due 
to stripe domains seen in Fig. [3] (b) . 

At 1.7 mT the splitting of the NSF reflectivities al- 
most vanishes. This qualitatively can be explained by a 
large angle between magnetization directions in neighbor- 
ing stripes. At the same time the SF reflectivity attains 
a maximum value manifesting a large projection of the 
mean magnetization component normal to the field as 
seen in Fig. [1] Hence, already qualitative analysis of the 
specular PNR and MOKE data let us conclude that the 
magnetization vectors in neighboring stripes in Fig. [3] (c) 
are not collinear with either the magnetic field or the EB 
axis and make quite a large angle between themselves. 
This angle can precisely be determined via the quantita- 
tive analysis of the complete scope of the PNR data. 

We have undertaken further PNR measurements at 
2.2 mT (not shown here) and 3.5 mT. At 3.5 mT in- 
creased splitting of the NSF and reduced SF reflectivities 
may indicate that the mean magnetization now turns to- 
wards the direction of the applied field. This conclusion, 
however, seems to contradict the MOKE data in Fig. [1] 
which do not show a substantial rotation of the mean 
magnetization within this interval of fields. Subsequent 
quantitative PNR analysis removes the contradiction be- 
tween the vector-MOKE observations and the PNR anal- 
ysis based on intuitive arguments. 

Further PNR measurements at 5.3 mT (not shown), 
6.7 mT and 7.1 mT (not shown) exhibit a continuous in- 
crease of the splitting of the NSF and a reduction of the 
SF reflectivities. The situations at 6.7 mT and 7.1 mT 
are comparable with the KM measurement in Fig. [3] (d) 
showing gradual shrinking and final collapse of the set of 
stripes with unfavorable magnetization. The SF reflec- 
tivities are already smaller as compared to the situation 
before the first reversal at 0.3 mT. Taking a closer look, 
one can also admit a number of particular details dis- 
tinguished in different plots for PNR, and in particular, 
those recorded at 0.3 and 7.1 mT. It is rather difficult 
to guess a physical meaning for most of the changes in 
the PNR q— dependencies. Nonetheless, the least square 
routine, as we shall see, allows us to infer a variation of 
a few field dependent parameters quantifying KM and 
MOKE observations. 

The major part of irrelevant parameters, e.g. those 
independent of the applied field are fixed via fitting the 
data collected at 361 mT, assuming the system at this 
field is in saturation. The maximum splitting seen in 
the last plot in Fig. [4] for the NSF reflectivities indicates 
that the magnetization is aligned along the applied field. 
Little SF reflectivity is observed due to a not perfect effi- 
ciency of the polarization device and taken into account 
in the subsequently applied least square routine. We also 
performed a measurement at a positive field of 3.1 mT 
in the backwards branch. It still shows a strong split- 
ting of NSF but an increased SF reflectivity compared to 
saturation indicating a tilted magnetization. 

The most complete information on the microscopic ar- 



0.04 
0.03 
0.02 
0.01 



exp. 



I++ 


/ 







f 



sim 



0.02 

Pi [A" 1 



II 







0.02, 


0.04 




Pi [A" 1 ] 






FIG. 5: (Color online) Experimental (left column) maps for 
the polarized neutron scattering intensity on a logarithmic 
scale, measured at a magnetic fields of 6.7 mT. The intensi- 
ties of the I + + , I non-spin- flip, and the intensities of the 
I + ~ , I~ + spin-flip cross-sections are plotted as a function of 
the angle of incidence at, and the scattering angle otf. Left 
column represents respective maps calculated in DWBA. 



rangement of magnetization is, however, obtained by an- 
alyzing not solely the specular reflection, but in accor- 
dance with off-specular scattering. Figure O displays ex- 
perimental data (left column) along with the results of 
theoretical simulations (right column) for all four scat- 
tering cross sections I + + , I , I + ~ , and I~ + collected 
into a set of maps. The scattering intensities are plot- 
ted as functions of the normal to the surface components 
Pi = k sin cti and pt = k sin a / of, correspondingly, in- 
coming, ki, and outgoing, kf, wave vectors impinging 
onto the surface at glancing angles of incidence, «j, and 
scattered at angles a/. The maps were obtained for all 
fields listed in Fig. [4] and in the text above but here we 
present only those constructed for one field at 6.7 mT 
and containing all features significant for the subsequent 
quantitative analysis. In the maps the specular refiec- 
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tion ridge runs along the diagonal, where Pi = Pf- At 
Pi 7^ Pf the scattering maps exhibit two other remarkable 
features. The first one is the intensity of Bragg diffrac- 
tion concentrated along curved lines « {2im/ h) cosx, 
where qu = (27r/A)(cos cti — cos a/), n denotes the or- 
der of diffraction, A is the period, and \ is the angle 
between the stripe axis and the normal to the reflec- 
tion planed. Bragg scattering occurs due to the periodic 
variation of the magnetization across the striped pattern 
and, in particular, due to the periodic alteration of the 
sign of (cos A7) coh in neighboring stripes. The diffracted 
intensity vanishes in saturation (not shown here) and 
should reach maximum values in the field range where 
antiparallel alignment of the magnetization in neighbor- 
ing microstripes is expected. The second feature is well- 
structured intensity of diffuse scattering observed at low 
angles of incidence oli and/or scattering a/. Both fea- 
tures are due to the lateral magnetization fluctuations 
on a scale smaller than the coherence range. 

In Fig. strong Bragg reflections for n = ±1 and 
weaker ones for n = ±2 can be recognized. The observa- 
tion of second-order Bragg reflections is quite a striking 
result. In the case of perfect alternation of magnetiza- 
tion projections in neighboring stripes of equal widths 
Bragg reflections of all even orders should be heavily 
suppressed due to the AF structure factor. Hence re- 
flection of the second order was never observed in our 
previous measurements^ carried out at % = 0. In the 
present case of x = 45° one of the KM images in Fig. 
e.g. at H w 7.0 mT, clearly indicates a difference in 
the widths of neighboring stripes. This difference vio- 
lates the cancellation law for the Bragg reflection of the 
even order—. Interestingly, second-order reflections can 
be observed not only when they are expected from the 
corresponding KM image in Fig. but rather at all fields 
along the ascending branch below saturation. In view 
of this, one should admit that the cancellation law re- 
quires a perfect reciprocity between magnetic moments 
of neighboring stripes. It can be violated not only be- 
cause of a difference in the stripe widths, but also due to 
non-perfect alternation of stripe magnetization projec- 
tions. This is particularly the case if the magnetic field 
is applied at an angle / with respect to the main 
symmetry axis. Then, in contrast to the symmetric case 
X = 0,— the external field tilts the magnetization vector 
by the angle A7 = fi\ in one set of stripes, or by A7 = /3 2 
in the other set. In the asymmetric case x there is 
no reason to expect that (3\ — —fo, while at /3i 7^ — (3% 
neither of the stripe magnetization projections perfectly 
alternate. 

Bragg reflections are observed in all four, SF and NSF, 
channels. The SF maps show a strong asymmetry with 
respect to the interchange pi with pf corresponding to 
parity between Bragg reflections with indexes n and — n. 
The asymmetry is explained by the birefringence^ in 
the mean optical potential and is accounted for within 
the framework of the distorted wave Born approximation 
(DWBA)22i2i. In the symmetric case (x = 0) /3i = -0 2 



and only the magnetization vector components collinear 
with the stripes alternate^. Then the SF diffraction is a 
result of the superposition of two effects: Bragg diffrac- 
tion due to periodical alternation of the scattering poten- 
tial and a homogeneous transverse magnetization which 
mixes up neutron spin states. This superposition is de- 
scribed in DWBA. In the present arrangement x / 
both in-plane projections of the stripe magnetization 
vectors alternate, providing either NSF and SF Bragg 
diffraction already in the Born approximation. However, 
an accurate balance between intensities in all channels 
as well as between specular and Bragg diffraction is only 
possible to account for accurately in DWBA. 

In Fig. Owe also observe diffuse off-specular scattering, 
which is due to random fluctuations A/3i and A/?2 around 
their mean values f3 1 and /3 2 , respectively The SF diffuse 
scattering is also strongly asymmetric and the asymme- 
try degree depends on the net magnetization projection 
onto the field guiding neutron polarization. In the / + ~~ 
map off-specular scattering intensity is mostly disposed 
at pj < pf, while in the I~ + one it is concentrated at 
Pi > pf. NSF diffuse scattering is, in contrast, symmet- 
ric. SF and NSF diffuse scattering together indicate that 
there are magnetization fluctuations of both, longitudinal 
and transverse components. Off-specular diffuse scatter- 
ing is strongly connected to the development of ripple 
domains and is most pronounced just below the first (at 
0.3 mT) and around the second coercive field (at 7.1 mT). 
At fields between the both coercive fields the diffuse SF 
intensity is much lower with minimum values at 3.5 mT, 
accounting for a much more regular domain state. Inter- 
estingly, in the descending branch at 3.1 mT we observe 
strongest diffuse SF scattering and no Bragg reflections. 
The reason for such behavior is nicely visualized in the 
KM measurement in Fig. |3fe) with strong ripple devel- 
opment and almost no stripe contrast accounting for a 
similar magnetization orientation in both stripe regions. 



IV. DATA ANALYSIS AND DISCUSSION 

As it has already been mentioned the PNR data, al- 
though containing a bulk of information, require a theo- 
retical model for their quantitative interpretation. Such 
a model founded on the vector-MOKE results and in par- 
ticular on the KM images which imply the existence of 
at least two types of hyper-domains comprising a num- 
ber of EB stripes. This means that the net magnetization 
vector 

M = wi(M 1 ) + W2(M 2 ) (5) 

is the sum of hyper-domain magnetization vectors (M±) 
and (M2) (see, Fig. \5§ weighted in accordance with the 
percentages w\ and W2 = 1 — toi of the sample area 
they cover. The directions of the vectors of the lo- 
cal magnetization M\p = (Afi.2) + AAfi.2 may vary 
as a function of the lateral coordinates x, y so that 
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FIG. 6: (Color online) Sketch of two hyper-domains separated 
by the domain wall (red dashed line) along one of the stripe 
boundary. The magnetization vectors of the hyper-domains 
are denoted by black arrows tilted by the angles 71,2 against 
the applied field. The red arrows indicate the local magne- 
tization vector tilted randomly and/or periodically by angles 
A7i,2 with respect to the hyper-domain magnetization direc- 
tions 



AAfi 2 describes local deviations of the local magneti- 
zation from their mean values (M\^ 2 ) averaged over each 
of the hyper-domains. Due to the strong anisotropy in 
the neutron coherency the averaging of the PNR signal, 
as it was pointed out above, should be performed in two 
steps. Firstly, the magnetization vectors are averaged 
over the coherence ellipsoids which are dramatically ex- 
tended along one axis but still are fitting into any of the 
hyper-domains. In the particular kinematics the long 
axis is parallel to the x— axis and the coherent averag- 
ing results in (Mi :2 (x : y)) co h = M 1:2 (y). The absolute 
values I Af 1,2(2/)] of these vectors determine the magnetic 
parts of the optical potentials and reflection amplitudes 
R 12 (y) for each type of hyper-domain. The amplitudes 
may still vary as a function of the y-coordinate and the 
equations for SF and NSF reflection intensities require 
secondly an additional incoherent averaging of the cor- 
responding cross sections along the y-direction. Those 
two types of averaging give access to not only the mean 
values (Mi, 2) = (M 1,2 (y))inc in Eq.(5) and, finally, to 
the net magnetization vector M, but also to the weights 
wi,2- Moreover, the least square routine provides us with 
quite a few parameters rather characterizing the domain 
model in great details. 

Our model assumes that the projections Mf 2 (a;, y) — 
M sat cos 7i )2 and Mf 2 (x,y) = M sat sin 71,2 of the mag- 
netization vectors are determined by the angles 71,2 = 
71,2(2;, y) = 71,2(2/) + A7i, 2 (a;,2/), where Aji, 2 (x,y) are 
describing deviations of the angles in directions of the 
vectors Mi j2 (x,y) from that of Mi_ 2 (y). The latter are 
tilted by 7^ 2 (y) against the a;— axis and are determined 
by the constrains (sin A7i ! 2) C oh = specific for each type 
of hyper-domains and the y-coordinate. The deviation in 



angles A7i,2, either random (ripple domains) and/or pe- 
riodic (stripes), reduce the absolute values |Mi j2 (y)| = 
M sat d. 2 (y) by the factors c h2 {y) = (cos(A7i i2 )) C oh < 
1. These factors generally depend on the y-coordinate. 
However, if the coherence length crosses a large num- 
ber of stripes and/or ripple domains this dependency 
is weak and can be neglected in the first approxima- 
tion so that only two parameters c\ and c 2 characteriz- 
ing reflection amplitudes are used in the fitting routine. 
Two other couples of parameters Ci 2 = (cos 7^ 2 (2/))inc 
and S 2 i :2 = (sin 2 7^ 2(2/))inc used to fit the data follow 
from the incoherent averaging of the PNR cross sections. 
Such an averaging accounts for fluctuations of the an- 
gles 7 y 1 2 (y) determined for different coherence ellipsoids. 
These fluctuations can be rather developed due to e.g. 
ripple domains which size is greater than the coherence 
length in the y-direction. 

The set of equations for NSF and SF reflectivities used 
in the data fitting are written as follow o 23 ' 24 : 

K{PuPf) = w 1 U 1 {P i ,P f ) + w 2 K 2 (P u Pf) 

Kx,2{Pi,p } ) = \{ [|i?+ 2 | 2 + |#i; 2 | 2 ][i + PiPf&iM 
+ i\RU 2 -\Ri.2n^ + Pf)Ci,2 

+2K(R+; 2 Rl 2 )P l P f S5 h2 }, (6) 

where P t = ±\P,\, P t = ±|P„-| with \Pi\ < 1 and |P;| < 1 
are efficiencies of the polarizer and analyzer, respectively. 
The complex reflection amplitudes R^ 2 — R(^ci, 2 ) are 
determined for two parameters c\ and c 2 . 

At each value of applied field the best fit of all 4 mea- 
sured reflection curves was obtained by varying 7 es- 
sential parameters: Ci,2, Ci ;2 , »5 2 i,2 = 1 — C 2 i,2 and 
w 2 = 1 — wi while keeping all others found from the fit 
at saturation where C1.2 = 1, Ci,2 = 1, S 2 i t2 = and 
u?i = 1. The quality of the fit is illustrated in Fig. [4] 
while the results are collected into Fig. [7j where the field 
variation of the parameters is presented. 

First of all one can admit that two types (Fig. [7Ja)) 
of hyper-domains, one with reduced (Fig. [7jb)) and the 
other with saturation magnetization, exist almost all 
over the range of the hysteresis loop. At low fields 
wi « w 2 w 0.5, i.e. the two domains states are equally 
populated. An increasing field H leads to a two step 
grow of the fraction wi (H) of hyper-domains with re- 
duced magnetization on the cost of those saturated until 
wi = 1 and w 2 = is reached just below the second co- 
ercive field. The reduction factor C\(H) of magnetization 
in the first type of hyper-domains plotted in Fig. [TJ^b) is 
mostly attributed to the angle between the magnetization 
vectors in neighboring stripes. It is not a monotonous 
function of field H and c\(H) has a minimum at H « 
2.2 mT where this angle has a maximum value. In ac- 
cordance to Fig. [JJc) the magnetization vectors of both 
types of hyper-domains rotate towards the magnetic field 
direction but with a different rate. The magnetization 
of "striped" hyper-domains approaches the applied field 
direction much faster than that of hyper-domains with 
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FIG. 7: (Color online) Results of fits to the PNR data with 
respect to the two hyper-domains and as function of the 
magnetic field, (a) Fraction Wi of hyper-domain 1, (b) re- 
duction factor ci = (cos(A7i)) coh , (c) parameters Ci,2 = 
(cos7 12 (y)) inc , (d) S 2 i, 2 = (sin 2 7 lj2 (t/))inc, and (e) disper- 
sion Ei, 2 = C 2 i,2 — C 12 . The solid symbols represent the 
results of fits to the PNR of the first hyper-domain and the 
hollow symbols that of the second hyper-domain. Lines in 
(b)-(d) represent respective results of the MOKE measure- 
ments for comparison. Please note that in (d) S 2 can only 
be achieved from PNR. Instead, the PNR data are compared 
to S as calculated from S from MOKE. In (e) the lines are 
guides to the eye. 



its magnetization close to saturation. Note that at the 
two lowest values H = 0.3 and H = 0.8 mT the mag- 
netization of both types of domains has a negative pro- 
jection onto the field direction. At the next measured 
field H = 1.7 mT the parameter C\ > 0, meaning that 
the magnetization vector ^-projection of the unsaturated 
hyper-domains now is positive. At the same field C2 < 0, 
i.e. the magnetization vectors of the saturated hyper- 
domains is tilted with respect to the field by angles of 
90° < 7 2 < 270°. Further increase of magnetic field pull 
both magnetization vectors to the field direction. 

Eq.(6) contains only S 2 i t 2 — (sin 2 7} plotted in 
Fig. 0(d), but not (sin 7) 2 . Therefore, as already men- 
tioned, with PNR alone we cannot determine the sense 
of rotation. The missing information is compensated due 



to the vector-MOKE measurements presented in Fig. [T] 
while PNR provides access to the dispersions Si ^ — 

C 2 i,2 — C 12) additional physical parameters plotted in 
Fig. E)(e). These quantities measure a degree of magneti- 
zation vector fluctuations in hyper-domains of the same 
type. Hence, £1 -C 1 all over the range of fields signifying 
on rather coherent rotation of the magnetization vector in 
different hyper-domains with reduced magnetization. For 
the other type of hyper-domains £2 , on the contrary, is 
rather high and reveals non-monotonous behavior. This 
means that the magnetization vector in different highly 
saturated hyper-domains is tilted by quite different an- 
gles 7 2 indicating the fact that there actually exist more 
types of hyper-domains with nearly saturated magnetiza- 
tion. Such hyper-domains can be distinguished on some 
of the KM images^!, but unfortunately, the quality of our 
present PNR data seems not allowing to introduce more 
parameters for their identification. 

However, our data, e.g. presented in Fig. [5] are well 
sufficient to infer the microscopic magnetic structure of 
hyper-domains and prove the model sketched in Fig. [5] 
Indeed, using parameters found from the fit of specular 
reflectivities we can identify the reference state of the 
system perturbed by periodic (stripe domains) and ran- 
dom (ripple domains) deviations of magnetization from 
its mean value in each type of hyper domains. Compar- 
ing intensity along the diffraction lines in the maps in 
Fig. [5] with specularly reflected intensity we determined 
the angles between magnetization directions in stripe do- 
mains and direction of magnetization of hyper-domain. 
For instance, at 6.7 mT it was found that f3 l = 30° 
and j3 2 — 0°, i.e. the angle between magnetization di- 
rections in neighboring stripes amounts 30°. The tilt 
angles was, actually, found by taking into account, that 
the intensity of the Bragg lines is reduced due to rip- 
ple domains. They cause magnetization fluctuations re- 
ducing the magnetic scattering contrast between stripe 
domains and simultaneously creating diffuse scattering 
seen in Fig. [5] The intensity of diffuse scattering allows 
us to quantify the amplitude of those fluctuations. In the 
case of H =6.7 mT ripple domains reduce stripe magne- 
tization for about 10% with respect to saturation. How- 
ever, fluctuations are not absolutely random but corre- 
lated over a few of neighboring stripes as is found from 
the extension of diffuse scattering in Fig. [5j The role 
of correlated magnetization fluctuations^ is not yet very 
clear, but they certainly strongly influence a formation of 
hyper-domains in the alternating EB systems and hence 
a scenario of the re-magnetization process. 



V. SUMMARY 

In summary, by using vector-MOKE magnetometry, 
Kerr microscopy and polarized neutron reflectomctry 
we studied the field induced evolution of the magne- 
tization distribution of a periodic pattern of alternat- 
ing exchange bias stripes when applying the magnetic 
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field at an angle of 45° in order to avoid the insta- 
bility with respect to the tilt of magnetization found 
for the easy axis configuration 13 . The data show that 
the re-magnetization process proceeds through different 
stages which are qualitatively described by comprehen- 
sive analysis of specular and off-specular PNR data. Be- 
side the formation of stripe-domains with alternated in- 
plane magnetization at small fields small ripple domains 
and two types large hyper-domains develop comprising a 
number of stripes. We developed a microscopic picture 



of the magnetic structure based on the results of our po- 
larized neutron studies which to some extent could be 
verified via Kerr microscopy. 
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given in Ref.— . 

This assumes that (cos A7) co h has the same value for dif- 
ferent coherence areas over the sample surface illuminated 
by the neutron beam. Alternatively, additional averaging 
of reflection intensity over values of the mean magnetiza- 
tion in different coherence spots has to be performed. 
At x = the stripes run perpendicular to the specular re- 
flection plane and qu « (2tto/A). It should be noticed that 
by rotating the sample by an angle \ 7^ Bragg reflections 
can be observe d 20 ' 21 at smaller angles until they merge to 
the reflection ridge at \ t/2. 

Higher order Bragg reflections are suppressed by the stripe 
form-factor and are not detectable in either the formes^ 3 -, 
or the present experiment. 



